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Abstract 

Phase-plane  techniques  are  used  to  analyze  a  quadratic  system  of  ordinary  differential 
equations  that  approximates  a  single  relaxation-time  system  of  partial  differential  equations 
used  to  model  transient  behavior  of  highly  elastic  non- Newtonian  liquids  in  shear  flow 
through  slit  dies.  The  latter  one-dimensional  model  is  derived  from  three-dimensional 
balance  laws  coupled  with  differential  constitutive  relations  well-known  by  rheologists. 
The  resulting  initial-boundary-value  problem  is  globally  well-posed  and  possesses  the  key 
feature:  the  steady  shear  stress  is  a  non-monotone  function  of  the  strain  rate.  Results  of 
the  global  analysis  of  the  quadratic  system  of  ode’s  lead  to  the  same  qualitative  features 
as  those  obtained  recently  by  numerical  simulation  of  the  governing  pde's  for  realistic 
data  for  polymer  melts  used  in  rheological  experiments.  The  analytical  results  provide 
an  explanation  of  the  experimentally  observed  phenomenon  called  spurt;  they  also  predict 
new  phenomena  discovered  in  the  numerical  simulation;  these  phenomena  should  also  be 
observable  in  experiments.  ■  • 
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1.  Introduction 


The  purpose  of  this  paper  is  to  analyze  novel  phenomena  in  dynamic  shearing  flows 
of  non-Newtonian  fluids  that  are  important  in  polymer  processing  [IT].  One  striking  phe¬ 
nomenon.  called  “spurt,”  was  apparently  first  observed  by  Vinogradov  et  a 1.  [19]  in  ex¬ 
periments  concerning  quasi  static  flow  of  monodispersive  polvisoprenes  through  capillaries 
or  equivalently  through  slit  dies.  They  found  that  the  volumetric  flow  rate  increased  dra¬ 
matically  at  a  critical  stress  that  was  independent  of  molecular  weight.  Until  recently, 
spurt  has  been  associated  with  the  failure  of  the  flowing  polymer  to  adhere  to  the  wall  [5]. 

The  focus  of  our  current  research  is  to  offer  an  alternate  explanation  of  spurt  and  related 
phenomena. 

Understanding  these  phenomena  has  proved  to  be  of  significant  physical,  mathemati¬ 
cal.  and  computational  interest.  In  our  recent  work  [12],  we  found  that  satisfactory  expla¬ 
nation  and  modeling  of  the  spurt  phenomenon  requires  studying  the  full  dynamics  of  the 
equations  of  motion  and  constitutive  equations.  The  common  and  key  feature  of  constitu¬ 
tive  models  that  exhibit  spurt  and  related  phenomena  is  a  non-monotonic  relation  between 
the  steady  shear  stress  and  strain  rate.  This  allows  jumps  in  the  steady  strain  rate  to  form 
when  the  driving  pressure  gradient  exceeds  a  critical  value;  such  jumps  correspond  to  the 
sudden  increase  in  volumetric  flow  rate  observed  in  the  experiments  of  Vinogradov  et  ai. 

The  governing  systems  used  to  model  such  one-dimensional  flows  are  analyzed  in  [12] 
by  numerical  techniques  and  simulation,  and  in  the  present  work  by  analytical  methods. 

The  systems  derive  from  fully  three-dimensional  differential  constitutive  relations  with  m- 
relaxation  times  (based  on  work  of  Johnson  and  Segalman  [S]  and  Oldroyd  [16]).  They 
are  evolutionary,  globally  well  posed  in  a  sense  described  below,  and  they  possess  discon¬ 
tinuous  steady  states  of  the  type  mentioned  above  that  lead  to  an  explanation  of  spurt. 

The  governing  systems  for  shear  flows  through  slit-dies  are  formulated  from  balance  laws 
in  Sec.  2. 

Specifically,  we  model  these  flows  by  decomposing  the  total  shear  stress  into  a  polymer 
contribution,  evolving  in  accordance  with  a  differential  constitutive  relation  with  a  single 
relaxation  time  and  a  Newtonian  viscosity  contribution  (see  system  ( JSO )  in  Sec.  2.).  The 
flows  can  aiso  be  modelled  by  a  system  based  on  a  differential  constitutive  law  with  two 
widely  spaced  relaxation  times  (see  system  ( JSO 2)  in  [13].)  but  no  Newtonian  viscosity 
contribution.  Numerical  simulation  [9.  12]  of  transient  flows  at  high  Weissenberg  (Debo¬ 
rah)  number  and  very  low  Reynolds  number  using  the  model  (JSO)  exhibited  spurt,  shape 
memory,  and  hysteresis:  furthermore,  it  predicted  other  effects,  such  as  latency,  normal 
stress  oscillations,  and  molecular  weight  dependence  of  hysteresis,  that  should  be  analysed 
further  and  tested  in  rheological  experiment. 

In  earlier  work.  Hunter  and  Slemrod  [7]  used  techniques  of  conservation  laws  to  study' 
the  qualitative  behavior  of  discontinuous  steady  states  in  a  simple  one-dimensional  vis-, 

coeiastic  model  of  rate  type  with  viscous  damping.  They  predicted  shape  memory  and  m _ L-J 

hysteresis  effects  related  to  spurt.  A  salient  feature  of  their  model  is  linear  instability  and - - 

loss  of  evoiutionaiity  in  a  certain  region  of  state  space.  _ 

The  objective  of  the  present  paper  is  to  develop  analytical  techniques,  the  results  of  _ 

which  verify  these  rather  dramatic  implications  of  numerical  simulation.  Based  on  scaling  Codes 
introduced  in  [12],  appropriate  for  the  highly  elastic  and  very'  viscous  polyisoprenes  used  in  ,nd/or 
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the  spurt-experiment,  we  are  led  to  study  the  following  pair  of  quadratic  autonomous  ordi¬ 
nary  differential  equations  that  approximates  the  governing  system  (J SO)  in  the  relevant 
range  of  physical  parameters  for  each  fixed  position  in  the  channel: 


(1.1) 


Here  the  doc  denotes  the  derivative  d/dt,  T  is  a  parameter  that  depends  on  the  driving 
pressure  gradient  as  well  as  position  x  in  the  channel,  and  £  >  0  is  a  ratio  of  viscosities. 
Systv..^  (1.1)  is  obtained  by  setting  a  =  0  in  the  momentum  equation  in  system  (JSO);  tins 
aDDroximation  is  reasonable  because  a  is  at  least  several  orders  of  magnitude  smaller  than 
e.  We  show  that  steady  states  of  system  ( JSO ),  some  of  which  axe  discontinuous  for  non¬ 
monotone  constitutive  relations,  correspond  to  to  critical  points  of  the  quadratic  system. 
We  deduce  the  local  characters  of  the  critical  points,  and  we  prove  that  system  (1.1)  has 
no  periodic  orbits  or  closed  separatrix  cycles.  Moreover,  this  system  is  endowed  with 
a  natural  Lyapunov-like  function  with  the  aid  of  which  we  are  able  to  determine  the 
global  dynamics  of  the  approximating  quadratic  system  completely  and  thus  identify  its 
globally  asymptotically  stable  crical  points  (i.e.  steady  states)  for  each  position  x.  This 
analysis  is  carried  out  in  Sec.  3  When  o,  the  ratio  of  Reynolds  to  Deborah  numbers,  is 
strictly  positive,  the  stability  of  discontinuous  steady  states  of  system  {JSO)  remains  to 
be  settled.  Recently,  Nohel,  Pego  and  Tzavaras  [15]  established  such  a  result  for  simple 
model  in  which  the  polymer  contribution  to  the  shear  stress  satisfies  a  single  differential 
constitutive  relation:  for  a  particular  choice,  their  model  and  system  {JSO)  with  a  >  0 
have  the  same  behavior  in  steady  shear.  Their  asymptotic  stability  result,  combined  with 
numerical  experiments  and  research  in  progress,  suggest  that  the  same  result  holds  for  the 
full  system  {JSO),  at  least  when  a  is  sufficiently  small. 

In  Sec.  4.. the  analysis  of  Sec.  3.  is  applied  to  eac^  noint  x  in  the  channel,  allowing 
us  to  explain  spurt,  shape  memory,  hysteresis,  and  other  eff  ects  originally  observed  in  the 
numerical  simulations  in  terms  of  a  continuum  of  phase  p  raits.  We  discuss  asymptotic 
expansions  of  solutions  of  systems  {JSO)  and  {JSOz)  of  Ref.  [13]  in  powers  of  c  that  enable 
us  to  explain  latency  (a  pseudo-steady  state  that  precedes  spurt).  The  asymptotic  analysis 
also  permits  a  more  quantitative  comparison  of  the  dynamics  of  the  two  models  when  c  is 
sufficient Iv  small.  In  Sec.  5.,  we  discuss  physical  implications  of  the  analysis,  particularly 
those  that  suggest  new  experiments.  In  Sec.  6..  we  draw  certain  conclusions.  Although  the 
anaivsis  in  this  paper  applies  only  to  the  special  constitutive  models  we  have  studied,  we 
expect  that  the  qualitative  features  of  our  results  appear  in  a  broad  class  of  non-Newtonian 
fluids.  Indeed,  numerical  simulation  by  Kolkka  and  Ierley  [10]  using  another  model  with  a 
single  relaxation  time  and  Newtonian  viscosity  exhibits  very  similar  character. 
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2.  A  Johnson-Segalman-Oldrovd  Model  for  Shear  Flow 

The  motion  of  a  fluid  under  incompressible  and  isothermal  conditions  is  governed  by 
the  balance  of  mass  and  linear  momentum.  The  response  characteristics  of  the  fluid  are 
embodied  in  the  constitutive  relation  for  the  stress.  For  viscoelastic  fluids  with  fading 
memory,  these  relations  specify  the  stress  as  a  functional  of  the  deformation  history  of  the 
fluid.  Many  sophisticated  constitutive  models  have  been  devised;  see  Ref.  [2]  for  a  survey. 
Oi  particular  interest  is  a  class  of  differential  models  with  m-relaxation  times,  derived  in  a 
three-dimensional  setting  in  Refs.  [12]  and  [13];  these  models  can  be  regarded  as  a  special 
cases  of  the  Johnson-Segalman  model  [Sjwnen  the  memory  function  is  a  linear  combina¬ 
tion  of  m-decaying  exponentials  with  positive  coefficients  or  of  the  Oldroyd  differential 
constitutive  equation  [16]. 

Essential  properties  of  constitutive  relations  are  exhibited  in  simple  planar  Pciseuille 
shear  flow.  We  study  shear  flow  of  a  non-Newtonian  fluid  between  parallel  plates,  located 
at  x  =  —h/2.  with  the  flow  aligned  along  the  y-axis,  symmetric  about  the  center  line,  and 
driven  by  a  constant  pressure  gradient  f .  We  restrict  attention  to  the  simplest  model  of  a 
single  relaxation-time  differential  model  that  possesses  steady  state  solutions  exhibiting  a 
non-monotone  relation  between  the  total  steady  shear  stress  and  strain  rate,  and  thereby 
reproduces  spurt  and  related  phenomena  discussed  below.  The  total  shear  stress  T  is 
decomposed  into  a  polymer  contribution  and  a  Newtonian  viscosity  contribution.  When 
restricted  to  one  space  dimension  the  initial-boundary  value  problem,  in  non-dimensional 
units  wi'h  distance  scaled  bv  h.  governing  the  flow  can  be  written  in  the  form  (see  Refs.  [9. 
12]): 

avt  -  <rr  =  £vxx  +  /  , 

<rt  —  (z  +  ik  =  -<x ,  (  jso ) 

Zt  4-  CTV z  =  — Z 

on  the  interval  [—1/2.0],  with  boundary  conditions 

u(  — 1/2,  t)  =  0  and  ur(0,f)  =  0  (BC) 

and  initial  conditions 

v(x.  0)  =  u0fr  ;  ,  cr(r.  0)  =  <7q(x)  ,  and  Z(x.  0)  =  Z0(x)  ,  on  —  1/2  <  x  <  0:  (JC) 

symmetry  of  the  flow  and  compatibility  with  the  boundary  conditions  requires  that 
vi)(  — 1/2)  =  0.  Uq(0)  =  0  and  a o(0)  =  0. 

The  evolution  of  <r .  the  polymer  contribution  to  the  shear  stress,  and  of  Z ,  a  quantity 
proportional  to  the  normal  stress  difference,  are  governed  by  the  second  and  third  equations 
in  system  (JSO).  As  a  result  of  scaling  motivated  by  numerical  simulation  and  introduced 
in  Ref.  [12],  there  are  only  three  essential  parameters:  a  is  a  ratio  of  Reynolds  number  to 
Deborah  number,  s  is  a  ratio  of  viscosities,  and  /  is  the  constant  pressure  gradient. 

When  £  =  0.  and  Z  +  1  >  0,  system  (JSO)  is  hyperbolic,  with  characteristics  speeds 
±((Z  4-  l)/aj1/2  and  0.  Moreover,  for  smooth  intial  data  in  the  hyperbolic  region  and 
compatible  with  the  boundary  conditions,  techniques  in  [IS]  can  be  used  to  establish 
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global  well-posedness  (in  terms  of  classical  solutions)  if  the  data  are  small,  and  finite¬ 
time  blow-up  of  classical  solutions  if  the  data  are  large.  If  e  >  0,  system  ( JSO )  for  any 
smooth  or  piece-wise  smooth  data;  indeed,  general  theory  developed  in  [15]  (see  Sec.  3 
and  particularly  Appendix  A)  yields  global  existence  of  classical  solutions  for  smooth 
initial  data  of  arbitrary  sire,  and  also  existence  of  almost  classical,  strong  solutions  with 
discontinuities  in  the  initial  velocity  gradient  and  in  stress  components;  the  latter  result 
allows  one  to  prescribe  discontinuous  initial  data  of  the  same  type  as  the  discontinuous 
steady  states  studied  in  this  paper. 

The  steady-state  solutions  of  system  (JSO)  play  an  important  role  in  our  discussion. 
Such  a  solution,  denoted  by  77,  a,  and  Z ,  can  be  described  as  follows.  The  stress  components 
a  and  Z  are  related  to  the  strain  rate  v z  through  the  relations 


(2.1) 


Therefore,  the  steady  total  shear  stress  T  :=  a  +  eux  is  given  by  T  =  w(vx),  where 


(2.2) 


The  properties  of  w,  the  steady-state  relation  between  shear  stress  and  shear  strain 
rate,  are  crucial  to  the  behavior  of  the  flow.  Bv  symmetry,  it  suffices  to  consider  s  >  0. 
For  ail  £  >  0.  the  function  w  has  inflection  points  at  s  =  0  and  s  =  a/3-  When  e  >  1/3, 
the  function  w  is  strictly  increasing,  but  when  £  <  1/8.  the  function  w  is  not  monotone. 
Lack  of  monotonicitv  is  the  fundamental  cause  of  the  non-Newtonian  behavior  studied  in 


this  paper:  hereafter  we  assume  that  £  <  1/S. 

The  graph  of  w  is  shown  in  Fig.  1.  Specifically,  w  has  a  maximum  at  s  =  s\/ 
and  a  minimum  at  s  =  sm,  where  it  takes  the  values  T \i  :=  w(s  V/)  and  Tm  w(sm) 
respectively.  As  £  —  1/S.  the  two  critical  points  coalesce  at  s  =  a/3". 

The  momentum  equation,  together  with  the  boundary  condition  at  the  centerline, 
implies  that  the  steady  total  shear  stress  satisfies  T  =  —  fx  for  every'  x  6  [—  0].  Therefore, 

the  steady  velocity  gradient  can  be  determined  as  a  function  of  x  by  solving 


vj(vx)  =  -fx  . 


(2.3) 


Equivalently,  a  steady  state  solution  vx  satisfies  the  cubic  equation  P(vx)  =  0,  where 

F(s)  :  =  £  s 3  -Ts2  +  (1  +  £)s-T  .  (2.4) 


The  steady  velocity  profile  in  Fig.  2  is  obtained  by  integrating  77r  and  using  the  boundary 
condition  at  the  wail.  However,  because  the  function  w  is  not  monotone,  there  might 
be  up  to  three  distinct  values  of  Ur  that  satisfy  Eq.  (2.3)  for  any'  particular  x  on  the 
interval  '  —  1/2.01.  Consequently.  77  x  can  suffer  jump  discontinuities,  resulting  in  kinks  in 
the  velocity  profile  (as  at  the  point  x.  in  Fig.  2).  Indeed,  a  steady  solution  must  contain 
such  a  jump  if  the  total  stress  T wa.u  =  f/2  at  the  wall  exceeds  the  total  stress  T m  at  the 
local  maximum  M  in  Fig.  1. 

Finally,  we  remark  that  the  flow  problem  discussed  here  can  also  be  modelled  by  a 
system  based  on  a  differential  constitutive  law  with  two  widely  spaced  relaxation  times 
but  no  Newtonian  viscosity  contribution  (see  system  (JSOn)  in  Sec.  2.  of  [13]);  with  an 
appropriate  choice  of  relevant  parameters,  the  resulting  proolem  exhibits  the  same  steady 
states  and  the  same  characteristics  as  (JSO). 
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Fig.  1:  Total  steady  shear  stress  T  vs.  shear  strain  rate  vx  for 
steady  flow.  The  case  of  three  critical  points  is  illustrated;  other 
possibilities  are  discusr  a  in  Sec.  3. 


3.  Phase  Plane  Analysis  for  System  ( JSO )  When  a  =  0 

When  a  is  not  zero,  numerical  simulation  developed  in  [9,  11,  12]  discovered  striking 
phenomena  in  shear  flow  and  suggested  the  analysis  that  follows.  A  great  deal  of  infor¬ 
mation  about  the  structure  of  solutions  of  system  (JSO)  can  be  garnered  by  studying 
a  quadratic  system  of  ordinary  differential  equations  that  approximates  it  in  a  certain 
parameter  range,  the  dynamics  of  which  is  determined  completely.  Motivation  for  this  ap¬ 
proximation  comes  from  the  following  observation:  in  experiments  of  Vinogradov  et  ai.  [19], 
a  is  of  the  order  10 ~ 1 2 ;  thus  the  term  avt  in  the  momentum  equation  of  system  (JSO) 
is  negligible  even  when  vt  is  moderately  large.  This  led  us  to  the  approximation  to  sys¬ 
tem  (JSO)  obtained  when  a  =  0. 

When  a  =  0,  the  momentum  equation  in  system  (JSO)  can  be  integrated  to  show 
that  the  total  shear  stress  T  :=  cr  +£VS  coincides  with  the  steady  value  T(x)  =  —  fx.  Thus 
T  =  T(x)  is  a  function  of  x  only,  even  though  cr  and  vc  are  functions  of  both  x  and  t.  The 
remaining  equations  of  system  (JSO)  yieid,  for  each  fixed  x,  the  autonomous,  quadratic, 
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Fig.  2:  Velocity  profile  for  steady  flow. 


p'ionar  system  of  ordinary  differential  equations 

'T-a 


<r  =  (Z  +  l)(: 


~  &  , 


(3.1) 


Here  the  dot  denotes  the  derivative  d/dt.  We  emphasize  that  for  each  /,  a  different 
dvnamicol  system  is  obtained  at  each  x  on  the  interval  [—1/2.0]  in  the  channel  because 
T  =  —fx.  Bv  svmmetrv,  we  may  focus  attention  on  the  case  T  >  0;  also  recall  from  Sec.  2 
that  s  <  l/S:  these  are  assumed  throughout.  The  dynamical  system  (3.1)  can  be  analyzed 
completely  bv  a  phase-plane  analysis  outlined  below;  the  reader  is  referred  to  Sec.  3  in  [13] 
for  further  details.  Here  we  state  the  main  results. 

The  critical  points  of  system  (3.1)  satisfy  the  algebraic  system 


(Z  +  l  +  e)  (|-l)+‘-  =  °. 

—  =  f=-:)~Z  =  0. 
e  T\T  J 


(3.2) 
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These  equations  define,  respectively,  a  hyperbola  and  a  parabola  in  the  a-Z  plane;  these 
curves  are  drawn  in  Fig.  3,  which  corresponds  to  the  most  comprehensive  case  of  three 
critical  points.  The  critical  points  are  intersections  of  these  curves.  In  particular,  critical 
points  lie  in  the  strip  0  <  a  <  T. 


Fig.  3:  The  pnase  plane  in  the  case  of  three  critical  points. 


Eliminating  Z  in  these  equations  shows  that  the  ^-coordinates  of  the  critical  points 
satisfy  the  cubic  equation  Q(cr/T)  =  0,  where 


Q{ 0  := 


-eu-i)  +  i+€ 


(£  ~  1)  +  £  . 


(3.3) 


A  straightforward  calculation  using  Eq.  (2.4)  shows  that 


P{V  r)  =  P 


'T -a 


=  --Q(a/T) 


(3.4) 


Thus  each  critical  point  of  the  system  (3.1)  defines  a  steady-state  solution  of  system  ( JSO ): 
such  a  solution  corresponds  to  a  point  on  the  steady  total-stress  curve  (see  Fig.  1)  at  which 
the  total  stress  is  T(x).  Consequently,  we  have: 


-8- 


Proposition  3.1: 

For  each  position  x  in  the  channel  and  for  each  s  >  0,  there  are  three  possibilities: 

(1)  there  is  a  single  critical  point  A  when  T  <  Tm; 

(2)  there  is  also  a  single  critical  point  C  if  T  >  T 

(3)  there  are  three  critical  points  A,  B ,  and  C  when  Tm  <  T  <  T m- 

For  simplicity,  we  ignore  the  degenerate  cases,  where  T  =  Tm  or  T  =  Tm,  in  which 
two  critical  points  coalesce. 

To  determine  the  qualitative  structure  of  the  dynamical  system  (3.1),  we  first  studv 
the  nature  of  the  critical  points.  The  behavior  of  orbits  near  a  critical  point  depends 
on  the  linearization  of  system  (3.1)  at  this  point,  i.e.,  on  the  eigenvalues  of  the  Jacobian 
matrix  J  associated  with  Eq.  (3.1),  evaluated  at  the  critical  point.  To  avoid  solving  the 
cubic  equation  Q(cr/T)  =  0,  the  character  of  the  eigenvalues  of  J  can  be  determined  from 
the  signs  of  the  trace  of  J  denoted  by  Tr  J,  the  determinant  of  J  denoted  by  Det  J.  and 
the  discriminant  of  J  denoted  by  Discrrn  J  at  the  critical  points.  We  omit  these  tedious 
calculations,  a  result  of  which  is  a  useful  fact:  at  a  critical  point,  e  Det  J  =  Q'(cr/T).  This 
relation  is  important  because  Q'  is  positive  at  .4  and  C  and  negative  at  B .  To  assist  the 
reader,  Fig.  3  shows  the  hyperbola  on  which  <7  =  0,  the  parabola  on  which  Z  =  0  [see 
Eqs.  (3.2)],  and  the  hyperbola  on  which  DiscrmJ  vanishes.  As  a  result  of  the  analysis 
above,  we  draw  the  following  conclusions: 

(1)  Tr  J  <  0  at  all  critical  points: 

(2)  Det  J  >  0  at  .4  and  C,  while  Det  J  <  0  at  f?;  and 

(3  i  DiscrmJ  >  0  at  .4  and  B ,  whereas  DiscrmJ  can  be  of  either  sign  at  C.  (For  tvpical 
values  of  c  and  T,  DiscrmJ  <  0  at  C;  in  particular.  DiscrmJ  <  0  if  C  is  the  only 
critical  point.  But  it  is  possible  for  DiscrmJ  to  be  positive  if  T  is  sufficientlv  close  to 

Tm.) 

Standard  theory  of  nonlinear  planar  dynamical  systems  (see.  e.g..  Ref.  [3.  Chap.  15] )  now 
establishes  the  local  characters  of  the  critical  points  A.  B ,  C  in  Proposition  3.1: 

Proposition  3.2: 

(1 )  .4  is  an  attracting  node  (called  the  classical  attractor); 

(2)  B  is  a  saddle  point: 

(3)  C  is  either  an  attracting  spiral  point  or  an  attracting  node  (called  the  spurt  attractor). 

The  next  task  is  to  determine  the  global  structure  of  the  orbits  of  system  (3.1).  In 
this  direction,  we  modify  an  argument  suggested  by  A.  Coppel  [4]  and  establish  the  cru¬ 
cial  result,  the  proof  of  which  involves  a  change  in  the  time  scale  and  an  application  of 
Bendixson’s  theorem: 

Proposition  3.3: 

System  (3.1)  has  neither  periodic  orbits  nor  separatrix  cycles. 

To  understand  the  global  qualitative  behavior  of  orbits,  we  construct  suitable  invariant 
sets.  In  this  regard,  a  crucial  tool  is  that  system  (3.1)  is  endowed  with  the  identity  (3.1) 

Jty+(Z  +  l)2}  =  -2  y  +  (Z  +  i)2  - 11  •  (3-5) 
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Tr.us  tne  function  \  (  Z  )  :—  a~  -r(Z  4-1)*  serves  as  a ^Lyapunov  function  for  the  dynamical 
system.  Notice  rha  identity  (3.5)  is  independent  of  T  and  s. 

Let  T  denote  the  circle  on  which  the  right  side  of  Eq.  (3.5)  vanishes,  and  let  Cr  denote 
the  circle  of  radius  r  centered  at  a  —  0  and  Z  =  — i.  i.e.  Cr  :=  {(cr,  Z)  :  V(cr,  Z)  =  r,r  >  0}; 
each  CP  is  a  level  set  of  V.  The  circles  T  and  C\  are  shown  in  Fig.  4,  which  corresponds 
*  -he  case  of  a  single  critical  point,  the  spiral  point  C .  Eq.  (3.5)  also  implies  the  critical 
cints  of  system  (3.1)  lie  on  T.  If  r  >  1,  T  lies  strictly  inside  Cr.  Consequently,  Eq.  (3.5) 
rows  that  the  dynamical  system  (3.1)  flows  inward  at  points  along  Cr.  Thus  the  interior 
f  Cr  is  a  positively  invariant  set  for  each  r  >  1.  Furthermore,  the  closed  disk  bounded  by 
.  which  is  the  intersection  of  these  sets,  is  also  positively  invariant.  Therefore  the  above 
rruneu:  establishes: 

roposition  3.4:  Each  closed  disk  hounded  by  the  circle  Cr,  r  >  1  is  a  positively 

■.variant  set  for  the  system  (3.1). 

The  above  results  combined  with  identification  of  suitable  invariant  sets  were  used  to 
nermire  the  global  structure  of  the  orbits  of  system  (3.1)  in  the  cases  of  one  and  three 
itical  points,  and  to  analyse  the  stable  and  unstable  manifolds  of  the  saddle  point  at  B. 
hese  results  are  shown  in  Figs.  5  and  6  and  summerized  in  the  following  result. 

roposition  3.5: 

he  5  asm  of  attraction  of  A.  i.e..  the  set  of  points  that  flow  toward  .4  as  t  — ►  oc.  comprises 
ti:  :se  points  on  the  same  side  of  the  stable  manifold  of  B  as  is  .4;  points  on  the  other  side 
a re  in  the  basin  of  attraction  of  C.  Moreover,  the  arc  of  the  circle  F  through  the  origin, 
between  5  and  its  reduction  B'  is  contained  in  the  basin  of  attraction  of  A.  In  particular, 
the  stable  manifold  for  B  cannot  cross  its  boundary  so  that  it  cannot  cross  T  between  B 
and  5'. 

.4.'.'  rub: mure  features  of  the  dynamics  of  system  (3.1)  (except  possibly  whether  C  is  a 
node  or  :j  focus  carry  over  to  one  that  approximates  the  system  ( JSO 2)  in  the  case  of  two 
widely  serarated  relaxation  times  (see  system  (4-3)  in  [13]). 


-10- 


z 


Fig.  4:  The  phase  plane  when  the  spurt  attractor  C  is  the  only 
critical  point. 

1.  Qualitative  Features  of  (J SO)  Based  on  Phase  Plane  Analysis 

The  discussion  that  follows  sketches  an  explanation  of  recent  numerical  simulations  of 
( JSO )  described  in  Refs.  [9,  12].  These  exhibited  several  effects  related  to  spurt:  latency, 
shape  memory,  and  hysteresis.  Fig.  7  shows  the  result  of  simulating  a  “quasi-static’’ 
loading  sequence  in  which  the  pressure  gradient  /  is  increased  in  small  steps,  allowing 
sufficient  time  between  steps  to  achieve  steady  flow  [9].  The  loading  sequence  is  followed 
by  a  similar  quasi-static  unloading  sequence,  in  which  the  driving  pressure  gradient  is 
decreased  in  steps.  The  initial  step  used  zero  initial  data,  and  succeeding  steps  used  the 
results  of  the  previous  step  as  initial  data.  The  resulting  hysteresis  loop  includes  the  shape 
memory  predicted  by  Hunter  and  Slemrod  [7j  for  a  simpler  model  by  a  different  approach. 
The  width  of  the  hysteresis  loop  at  the  bottom  can  be  related  directly  to  the  molecular 
weight  of  the  sample  [9]. 

We  explain  spurt,  shape  memory,  hysteresis  and  latency.  We  consider  experiments 
of  the  following  type:  the  flow  is  initially  in  a  steady  state  corresponding  to  a  forcing 
/0,  and  the  forcing  is  suddenly  changed  to  /  =  /„-(-  A/.  We  call^  this  process  “loading’’ 
(resp.  “unloading")  if  A/  has  the  same  (resp.  opposite)  sign  as  /0.  The  initial  flow  can 
be  described  by  specifying,  for  each  channel  position  r,  whether  the  flow  is  at  a  classical 
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Fig.  5:  The  orbit  through  origin  when  the  spurt  attractor  C  is  the 
only  critical  point. 

attractor  A  (x  is  a  “classical  p  'nt”)  or  a  spurt  attractor  C  (x  is  a  “spurt  point”)  for  the 
system  (3.1)  with  T  =  -fQx.  We  shall  say  that  any  point  lying  on  the  same  side  of  the 
stable  manifold  of  B  as  is  A  lies  on  the  “classical  side”;  points  lying  on  the  other  side  are 
said  to  be  on  the  “spurt  side.”  The  outcome  of  the  experiment  depends  on  the  character  of 
the  phase  portrait  with  T  =  -fx.  To  determine  this  outcome,  we  need  only  decide  when 
a  classical  point  becomes  a  spurt  point  or  vice  versa. 

The  principle  mathematical  properties  of  the  dynamical  system  (3.1)  that  determine 
the  outcome  of  loading  and  unloading  experiments  are  embodied  in  the  following  conse¬ 
quence  of  the  phase  plane  analysis. 

Proposition  4.1:  _ 

(1)  A  classical  point  .i0  for  the  initial  forcing  fQ  lies  in  the_domain  of  attraction  of  the 
classical  attractor  A  r  f,  provided  that  A  exists  (i.e.,  \fx\  <  T M); 

(2)  A  spurt  point  C  o  for  ,ne  initial  forcing  f0  lies  in  the  domain_of  attraction  of  the  spurt 
attractor  C  for  /  unless  (a)  C  does  not  exist  (i.e.,  \fx\  <  TJ;  or  (b)  C  lies  on  the 
classical  side  of  the  stable  manifold  of  the  saddle  point  B  for  f . 

Consider  starting  with  /0  =  0  and  loading  to  J  >  0.  Thus  the  initial  state  for  each  x 
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Fig.  6:  Phase  portrait  in  the  case  of  three  critical  points,  with  C 
being  a  spiral. 

lies  at  the  origin  <7  =  0,  Z  =  0.  Then  according  to  4.1(1)  above,  each  x  €  [-1/2,0]  such 
that  ~f\x\  <  T m  is  a  classical  point,  while  the  x  for  which  /|x|  >  TM  are  spurt  points 
(because  there 'is  no  classical  attractor).  Consequently,  we  draw  two  conclusions: 

Proposition  4.2:  _  _  _ 

(a)  If  die  forcing  is  subcritical  (i.e.,  f  <  /crit  :=  2 TM),  the  asymptotic  steady  flow  is 

entirely  classical.  _  _ 

(b)  If  the  forcing  is  supercritical  (f_>  fern),  there  is  a  single  kink  in  the  velocity  profile 
(see  Fig.  2),  located  at  x.  =  -T  m  /  f ;  those  x  €  [-1/2. x.),  near  the  wall  are  spurt 
points,  whereas  x  €  (x.,0],  near  the  centerline,  axe  classical. 

The  solution  in  case  (b)  can  be  described  as  “top  jumping'’  because  the  stress  T,  =  TlV* 
at  the  kink  is  as  large  as  possible,  and  the  the  kink  is  located  as  close  as  possible  to  the 

wall.  _  _  _  . 

Next,  consider  increasingjhe  load  from  /„  >  0  to  /  >  A  point  x  that  is  classical 

for  /0  remains  classical  for  /  unless  there  is  no  classical  attractor  for  T  -  -fx.  i.e., 
J\i\  >  Tv/  A  spurt  point  x  for  /0,  on  the  other  hand,  is  always  a  spurt  point  for  /.  As 
a  result,  a  point  in  x  in  the  channel  can  change  only  from  a  classical  attractor  to  a  spurt 
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Fig.  7:  Hysteresis  under  cyclic  load:  normalized  throughput  6 Q 
vs.  wall  shear  stress  Twail  [9]- 

attractor,  and  then  only  if  f\x\  exceeds  T,vf.  When  /  is  chosen  to  be  supercritical,  loading 
causes  the  position  x,  of  the  kink  in  Fig.  2  to  move  away  from  the  wall,  but  only  to  the 
extent  that  it  must:  a  single  jump  in  strain  rate  occurs  at  xm  =  — T,vf//,  where  the  total 
stress  is  T .  =  2\vf-  These  conclusions  axe  valid,  in  particular,  for  a  quasi-static  process  of 
gradually  increasing  the  load  from  fQ  =  0  to  f  >  /crit. 

Now  consider  unloading  from  /0  >  0  to  /  <  /0;  assume,  for  the  moment,  that  /  is 
positive.  Here,  the  initial  steady  solution  need  not  correspond  to  top  jumping.  For  this 
type  of  unloading,  a  point  x  that  is  classical  for  f0  always  remains  classical  for  /:  the 
classical  attractor  for  /  exists  because  / |x|  <  /0 jxj.  By  contrast,  a_spurt  point  x_for  /0 
can  become  classical  at  /.  This  occurs  if:  (a)  the  total  stress  T  =  —fx  falls  below  Tm;  or 
(b)  the  spurt  attractor  C0  for  T  =  —f0x  lies  on  the  classical  side  of  the  stable  manifold  of 
the  saddle  point  B  for  T  =  —  fx  (see  Proposition  4.1(2b)). 

Combining  the  analysis  of  loading  and  unloading  leads  to  the  following  summary  of 
quasi-static  cycles  and  the  resulting  flow  hysteresis. 

Kinks  move  away  from  the  wall  under  top  jumping  loading ;  they  move  toward  the  wall 
under  bottom  jumping  unloading;  otherwise  they  remain  fixed.  The  hysteresis  loop  opens 
from  the  point  at  which  unloading  commences;  no  part  of  the  unloading  path  retraces  the 
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loading  path  until  point  d.  in  Fig.  7. 

To  explain  the  latency  effect  that  occurs  during  loading,  assume  that  e  is  small.  It  is 
readily  seen  thac  the  total  stress  T m  at  the  the  local  maximum  M  is  1/2  +  0(e),  while 
the  local  minimum  m  corresponds  _to  a  total  stress  Tm  of  2x/£[l  4-  0(e)].  Furthermore, 
for  x  such  that  T(x)  =  0(1),  a  =  T  -f-  0(e)  at  an  attracting  node  at  A ,  while  a  =■  0(e) 
at  a  spurt  attractor  C  (which  is  a  spiral).  Consider  a  point  along  the  channel  for  which 
T[x)  >  T.v/.  so  that  the  only  critical  point  of  the  system  (3.1)  is  C,  and  suppose  that  that 
T  <  1.  Then  the  evolution  of  the  system  exhibits  three  distinct  phases,  as  indicated  in 
Fig.  6:  an  initial  “Newtonian*’  phase  (O  to  N );  an  intermediate  “latency”  phase  (jV  to  5); 
and  a  final  “spurt"’  phase  (5  to  O). 

The  Newtonian  phase  occurs  on  a  time  scale  of  order  £,  during  which  the  svstem 
approximately  follows  an  arc  of  a  circle  centered  at  cr  =  0  and  Z  =  —1.  Having  assumed 
that  T  <  1,  Z  approaches 

ZN  =  (1  -  T  )  j  -  1  (4.1) 

as  a  rises_to  the  value  T.  (If,  on  the  other  hand,  T  >  1,'the  circular  arc  does  not  extend 
as  far  as  T,  and  cr  never  attains  the  value  T;  rather,  the  system  slowly  spirals  toward  the 
spurt  attractor.  Thus  the  dynamical  behavior  does  not  exhibit  distinct  phases.) 

The  latency  phase  is  characterized  by  having  a  =  T-rO(e),  so  that  a  is  nearly  constant 
and  Z  evolves  approximately  according  to  the  differential  equation 


i  = 


T 

Z  +  l 


-Z  . 


(4.2) 


Therefore,  the  shear  stress  and  velocity  profiles  closely  resemble  those  for  a  steady  solution 
with  no  spurt,  but  the  solution  is  not  truly  steady  because  the  normal  stress  difference 
Z  sr.il  changes.  Integrating  Eq.  (4.2)  from  Z  =  Zy  to  Z  —  —1  determines  the  latency 
period.  This  period  becomes  indefinitely  long  when  the  forcing  decreases  to  its  critical 
value:  thus  the  persistence  of  the  near-steady  solution  with  no  spurt  can  be  very  dramatic. 
The  solution'  remains  longest  near  point  L  where  Z  =  —1  +  T.  This  point  may  be  regarded 
as  the  remnant  of  the  attracting  node  .4  and  the  saddle  point  B.  Eventually  the  solution 
enters  the  spurt  phase  and  tends  to  the  critical  point  C.  Because  C  is  an  attracting  spiral, 
the  stress  oscillates  between  the  shear  and  normal  components  while  it  approaches  the 
steady  state. 

Asymptotic  analysis  carried  out  in  Sec.  6  of  [13]  shows  that  when  s  is  sufficiently 
small,  system  ( JSOz )  of  [13]  has  the  same  asymptotic  properties  as  system  ( JSO ).  Thus 
system  (  JSO)  approximates  ( JSOi )  quantitatively  as  well  as  qualitatively. 
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5.  Physical  Implications 

One  of  the  widely  accepted  explanations  of  spurt  and  similar  observations  is  that  the 
presence  of  the  wall  affects  the  dynamics  of  the  polymer  system  near  the  wall.  Conceivably, 
there  could  be  a  variety  of  “wall  effects,”  the  most  obvious  is  the  loss  of  chemical  bond 
between  wall  and  fluid,  or  wall  slip  [5],  Perhaps  the  most  distinguishing  feature  of  our 
alternative  approach  is:  it  predicts  that  spurt  stems  from  a  material  property  of  the 
polymer  and  is  not  related  to  any  external  interaction.  The  spurt  layer  forms  at  the  wall 
in  situations  such  as  top  jumping  because  the  stresses  are  higher  there;  for  the  same  reason, 
of  course,  is  chemical  bonds  would  break  at  the  walljhowever,  our  approach  predicts  that 
the  layer  of  spurt  points  spreads  into  the  interior  of  the  channel  on  continued  loading 
Layer  thickness  is  predicted  to  grow  continuously  in  loading  to  a  thickness  that  should  be 
observable,  provided  secondary  (two-dimensional)  instabilities  do  not  develop. 

Our  analysis  suggests  other  ways  in  which  experiments  might  be  devised  to  verify 
the  dependence  of  spurt  on  material  properties:  (i)  produce  multiple  kinks  with  spurt 
layer  separated  from  the  wall,  (ii)  produce  hysteresis  in  flow  reversal  (Fig.  9).  Our  model 
predicts  circumstances  under  which  a  different  path  can  be  followed  in  sudden  reversal  of 
the  flow  than  would  be  followed  by  a  sequence  of  solutions  in  which  the  pressure  gradient 
is  reduced  to  zero  and  reloaded  again  (with  the  opposite  sign)  to  a  value  of  somewhat 
smaller  magnitude.  Such  behavior  does  not  seem  likely  to  be  explainable  by  a  wall  effect. 

The  most  important  and  perhaps  the  easiest  experiment  to  perform  to  verify  our  the¬ 
ory  is  to  produce  latency.  Our  analysis  predicts  long  latency  times  for  data  corresponding 
to  realistic  material  data:  no  sophisticated  timing  device  would  be  required,  nor  would  the 
onset  of  the  instability  be  hard  to  identify.  The  increase  in  throughput  is  predicted  to  be 
so  dramatic  that  simple  visual  inspection  of  the  exit  flow  would  probably  be  sufficient. 

6.  Conclusions 

Although  our  analysis  applies  only  to  the  special  constitutive  models  we  have  studied, 
we  expect  that  the  qualitative  features  of  our  results  appear  in  a  broad  class  of  non- 
Newtonian  fluids.  Our  analysis  has  identified  certain  universal  mathematical  features  in 
the  shear  flow  of  viscoelastic  fluids  described  by  differential  constitutive  relations  that 
give  rise  to  spurt  and  related  phenomena.  The  key  feature  is  that  there  axe  three  widely 
separated  time  scales,  each  associated  with  an  important  non-dimensional  number  (a.  £, 
and  1.  respectively),  when  scaled  by  the  dominant  relaxation  time.  A-1.  Each  of  these 
time  scales  can  be  associated  with  a  particular  equation  in  system  ( JSO )  [13].  The  key 
to  understanding  the  dynamics  of  such  systems  is  fixing  the  location  of  the  discontinuity 
in  the  strain  rate  induced  by  the  non-monotone  character  of  the  steady  shear  stress  vs. 
strain  rate. 
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